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Effective field theories (EFTs) have been widely used as a framework in order to place constraints 
on the Planck suppressed Lorentz violations predicted by various models of quantum gravity. There 
are however technical problems in the EFT framework when it comes to ensuring that small Lorentz 
violations remain small — this is the essence of the “naturalness” problem. Herein we present 
an “emergent” space-time model, based on the “analogue gravity” programme, by investigating a 
specific condensed-matter system that is in principle capable of simulating the salient features of an 
EFT framework with Lorentz violations. Specifically, we consider the class of two-component BECs 
subject to laser-induced transitions between the components, and we show that this model is an 
example for Lorentz invariance violation due to ultraviolet physics. Furthermore our model explicitly 
avoids the “naturalness problem”, and makes specific suggestions regarding how to construct a 
physically reasonable quantum gravity phenomenology. 


Introduction: The purpose of quantum gravity phe¬ 
nomenology (QGP) is to analyze the physical conse¬ 
quences arising from various models of quantum gravity 
(QG). One hope for obtaining an experimental grasp on 
QG is the generic prediction arising in many (but not 
all) quantum gravity models that ultraviolet physics at 
or near the Planck scale, Mpi anc k = 1.2 x 10 19 GeV/c 2 , 
(or in some models the string scale), typically induces vi¬ 
olations of Lorentz invariance (LI) at lower scales 13,0 
Interestingly most investigations, even if they arise from 
quite different fundamental physics, seem to converge 
on the prediction that the breakdown of Lorentz invari¬ 
ance (LI) can generically become manifest in the form 
of modified dispersion relations exhibiting extra energy- 
dependent or momentum-dependent terms, apart from 
the usual quadratic one occurring in the Lorentz invari¬ 
ant dispersion relation. In particular one most often con¬ 
siders Lorentz invariance violations (LIV) in the boost 
subgroup, leading to an expansion of the dispersion rela¬ 
tion in momentum-dependent terms, 

w" = Wq + (1 + 772) c 2 k 2 + 774 ^ Mliv ) k 4 + ... , (1) 

where the coefficients if n are dimensionless (and possibly 
dependent on the particle species considered), and we 
have restricted our expansion to CPT invariant terms 
(otherwise one would also get odd powers in k). The 
particular inertial frame for these dispersion relations is 
generally specified to be the frame set by Cosmological 
Microwave Background (CMB), and Mliv is the scale of 
Lorentz symmetry breaking which furthermore is gener¬ 
ally assumed to be of the order of -Mpi anc k. 

Although several alternative scenarios have been con¬ 
sidered in the literature in order to justify modified kine¬ 
matics of the kind of Eq. o, so far the most commonly 


explored avenue is an effective field theory (EFT) ap¬ 
proach. In the present article we wish to focus on the 
class of non-renormalizable EFTs with Lorentz violations 
associated to dispersion relations like Eq. © Relaxing 
our CPT invariance condition this class would include 
the model developed in §j], and subsequently studied by 
several authors, where an extension of QED including 
only mass dimension five Lorentz-violating operators was 
considered. (That ansatz leads to order k 3 LI and CPT 
violating terms in the dispersion relation.) Very accu¬ 
rate constraints have been obtained for this model using 
a combination of experiments and observations (mainly 
in high energy astrophysics). See e.g. 0]. 

In spite of the remarkable success of this framework 
as a “test theory”, it is interesting to note that there 
are still significant open issues concerning its theoreti¬ 
cal foundations. This is often referred to as the nat¬ 
uralness problem and can be expressed in the follow¬ 
ing way. Looking back at our ansatz m we can see 
that the lowest-order correction, proportional to 772, is 
not explicitly Planck suppressed. This implies that such 
term would always be dominant with respect to the 
higher order ones and grossly incompatible with obser¬ 
vations (given that we have very good constraints on 
the universality of the speed of light for different ele¬ 
mentary particles). Following the observational leads it 
has been therefore often assumed either that some sym¬ 
metry (other than Lorentz invariance) enforces the 772 
coefficients to be exactly zero, or that the presence of 
some other characteristic EFT mass scale ft, -C Mpi anc k 
(e.g., some particle physics mass scale) associated with 
the Lorentz symmetry breaking might enter in the low¬ 
est order dimensionless coefficient 772 — which will be 
then generically suppressed by appropriate ratios of this 
characteristic mass to the Planck mass: 772 oc ( ft/Mp\) a 
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where er > 1 is some positive power (often taken as one or 
two). If this is the case then one has two distinct regimes: 
For low momenta p/(Mpjc) -C {p/Mp\) a the lower-order 
(quadratic in the momentum) deviations in Q will dom¬ 
inate over the higher-order ones, while at high energies 
p/(Mp\c ) ( ii/Mp\Y the higher order terms will be 

dominant. 

The naturalness problem arises because such a scenario 
is not well justified within an EFT framework; in other 
words there is no natural suppression of the low-order 
modifications in these models. In fact we implicitly as¬ 
sumed that there are no extra Planck suppressions hidden 
in the dimensionless coefficients r] n with n > 2. EFT can¬ 
not justify why only the dimensionless coefficients of the 
n < 2 terms should be suppressed by powers of the small 
ratio n/Mp\. Even worse, renormalization group argu¬ 
ments seem to imply that a similar mass ratio, fi/Mp\ 
would implicitly be present also in all the dimensionless 
n > 2 coefficients - hence suppressing them even fur¬ 
ther, to the point of complete undetectability. Further¬ 
more it is easy to show 4] that, without some protecting 
symmetry, it is generic that radiative corrections due to 
particle interactions in an EFT with only Lorentz vio¬ 
lations of order n > 2 in o for the free particles, will 
generate n = 2 Lorentz violating terms in the dispersion 
relation^ which will then be dominant. Observational ev¬ 
idence [lj, |5j suggests that for a variety of standard model 
particles I 772 1 < 10 -21 . Naturalness in EFT would then 
imply that the higher order terms are at least as sup¬ 
pressed as this, and hence beyond observational reach. 

In order to contribute to this debate, we have cho¬ 
sen a rather unconventional path: We have investigated 
a condensed matter analog model (AM) of an emergent 
spacetime Jg), that reproduces the salient features of the 
the non-renormalizable EFT with LIV adopted in quan¬ 
tum gravity phenomenology studies. In particular we 
looked for a condensed matter system characterized by 1) 
relativistic kinematics for the low-energy quasi-particles 
2) presence of Lorentz violation at high energies due to 
the underling microscopic structure 3) coexistence on the 
same background of more than one quasi-particle species 
(in order to be able to see the presence of deviations from 
LI at order k 2 ). 

Standard Bose-Einstein condensates (BEC) are in this 
sense interesting systems as they fulfill the first two of the 
above requirements (see e.g. |7|, [s| ): their excitations can 
be described at low energy as relativistic phonons prop¬ 
agating on a geometrical background, and at higher or¬ 
der their dispersion relations show modifications of order 
k A (this is the so called Bogoliubov dispersion relation, 
which has been experimentally confirmed 0 )- Unfortu¬ 
nately in a single BEC system there is only one species of 
phonon, and hence it is impossible to address the ques¬ 
tion of naturalness. For this reason we chose to investi¬ 
gate the energy dependent behavior of quasi-particles in 
a 2-component Bose-Einstein condensate. For a specific 


choice of parameters such a system allows both a mass¬ 
less and a massive quasi-particle (see 0 ) , which share 
the same relativistic causal structure in the low energy 
limit. It is then natural to investigate the expected vio¬ 
lation of Lorentz invariance as the high energy regime of 
the theory is probed. (For a detailed discussion, see mt) 
Sound waves in 2-component BECs: The basis for our 
model is an ultra-cold dilute atomic gas of N bosons 
in two single-particle states |A) and | B). For example 
one could consider a two-component condensate of 87 Rb 
atoms in different hyperfine levels (see e.g. 00 The 
two states have slightly different energies, which permits 
us, from a theoretical point of view, to keep m.A 7 ^ mg, 
even if in experimentally realizable situations mA ~ ms- 
There are three atom-atom coupling constants, U aa-, 
Ubb , and Uab , and for our purposes it is essential to 
include an additional coupling A that drives transitions 
between the two single-particle states. For temperatures 
at or below the critical BEC temperature, almost all 
atoms occupy the respective ground states |A), \B) and it 
is hence meaningful to adopt the mean-field description 
for these modes. Ignoring back reaction effects of the 
quantum fluctuations one then obtains a pair of coupled 
Gross-Pitaevskii equations (GPE) 


iTidt'S’i 


' r -,2 


V + Vi — p,i 


2 mi 

+u u i^r+c/yi^-. 


+ A , (2) 


where ( i,j ) —> ( A,B ) or (i,j) —> (B,A) and is the 
classical wave function of the condensate (US'). Now con¬ 
sider small perturbations (sound waves) in the conden¬ 
sate cloud. The excitation spectrum is obtained by lin¬ 
earizing around some background, and after a straight¬ 
forward analysis leads to a differential equation for the 
(rescaled) perturbations in the phases, 6 = E~ x ^ 2 8 = 
S -1 / 2 [6Ui, $bi] T , where E is a 2 x 2 matrix constructed 
from the atomic couplings Uij 0 

Hydrodynamic limit: If one ignores the effect of the 
quantum potential then one obtains 


d 2 0 = -dt (i tfo • Vfl) - V ■ (Vo 1 0 ) 


+V' 


(eg - v 0 1 v 0 ) w +n 2 e, ( 3 ) 


where Cg and fl 2 are 2x2 symmetric matrices con¬ 
structed from the parameters appearing in the GPE, and 
vq is the common flow velocity of the background con¬ 
densates. If [Cg, U 2 ] = 0 then decomposition onto the 
eigenstates of the system results in a pair of independent 
“curved spacetime” Klein -Gordon equations 


\T~9 \/n 


d a ffi/n (si/ii) ab 9b6j / n j + Wj/n 6y n — 0, 

(4) 
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where the “acoustic metrics” are given by 


(Sl/Il)ab oc 


r (2 .,2^ 

1 -* T 1 

1 -vo 

- ( C I/II - v o) 

T 0 

r 5> 

1 

1 

| I dxd 


(5) 


and where the overall conformal factor depends on the 
spatial dimension d. The metric components depend only 
on the background velocity iTo, the background densities 
poi, and the speeds of sound c 2 / r] for the two eigenmodes. 
These are given by the eigenvalues of the matrix C 


„ 2 _ tr[Cg] ± \/tr [Cq] 2 — 4det[Cg] 

e i/n —- 2 -' l 0 ' 

The speed of sound in the AM takes on the role of the 
speed of light. The matrix Q 2 can be shown to have 
zero determinant, and so the eigenfrequencies of the two 
phonon modes are: uif = 0 and c/ 2 = tr[0 2 ]. The masses 
of the modes are then defined as m^ n = 
and thus the AM corresponds to one massless particle 
mi = 0 and one massive particle. They both “experi¬ 
ence” the same space-time if the sound speeds are equal, 
which requires tr[Gg] 2 — 4det[Gg] = 0, i.e. cj = Cu = 
cq = tr[Gg]/2. Hence we now have an AM representing 
both massive and massless particles, propagating on the 
same background at low energies. Let us now extend the 
analysis to high-energies and explore the structure of the 
corresponding LIV. 

QGP beyond the hydrodynamic limit: Starting from 
the GPE, we now linearize around a uniform condensate 
and set the background velocity to zero, vq = 0, but 
retain the quantum pressure term. The equationjor the 
phase perturbations in momentum space is now 0 


u> 2 9 


jV“+X k 2 [D k 2 + A] a/S + X fc 2 j 9 
H(k 2 ) 9 , (7) 


where S, X, D , and A are additional 2x2 symmetric ma¬ 
trices constructed from the parameters appearing in the 
GPE. The perturbation spectrum obeys the generalized 
Fresnel equation: 


in k , as the dispersion relation is invariant under CPT. 
This is by no means a surprising result, because the un¬ 
derlying GPE 0 is also invariant under CPT. 

It is now useful to define the symmetric matrices 
C 2 = C% + AC 2 and AG 2 = X^AX 1 / 2 which de¬ 
scribe how the speed of sound is modified, and also to 
define V 2 = 2X 1 / 2 S“ 1 X 1 / 2 and X 2 = 2 X 1 I 2 DX 1 ' 2 . 
All three of AG 2 , X 2 , and Z 2 are explicitly suppressed 
by powers of the mass of the fundamental constituents 
that condense to form the BEC. Note that all the rele¬ 
vant matrices have been carefully symmetrized, and note 
the important distinction between Gq and G 2 . Now de¬ 
fine c 2 = ^tr[G 2 ], which approaches the speed of sound 
c 2 —> Cq, in the hydrodynamic limit. The second and 
fourth order coefficients in the dispersion relations m 
are jlH ]: 


d^i/n 

r 2 

( 2tr[fl 2 Gg] — tr [ST 2 ] tr[G 2 ] 1 

dfc 2 

k->0 

± 1 tr[G 2 ]tr[fi 2 ] J 


dV/n 


d(fc 2 


c 2 (l ± 772); 
1 


( 11 ) 


k^r 0 


2 

±2 


tr [Z 2 ] ± tr [Z 2 ] 

tr[fl 2 Gg] — tr[fl 2 ] tr[Gg] 
tr[fl 2 ] 


tr[X 2 


± tr[G 2 ] 2 -jtdet[G 0 2 ] ^ tr[G 2 ] 2 ^ 2 


tr[fl 2 ] 
= 2?74 (^/Mliv)‘ 


tr[fl 2 


( 12 ) 


Lorentz violations from UV physics: In order to obtain 
LIV purely due to ultraviolet physics, we demand exact 
Lorentz invariance in the hydrodynamic limit. In other 
words, we require all terms in the equations GH> and m 
which might otherwise survive in the hydrodynamic limit 
to be set to zero. The constraints we obtain are: 

G1 : tr[G 2 ] 2 - 4det[G 2 ] = 0; (13) 

G2 : 2tr[f2 2 Gg] - tr[0 2 ]tr[G^] = 0 . (14) 


det{w 2 I — H(k 2 )} = 0, (8) 


and the dispersion relations for the phonon modes are 


tr[U(fc 2 )] ± ,/tr[Lf(fc 2 )] 2 - 4 det [H(k 2 )} 


- 1/11 2 

A Taylor-series expansion gives 


(9) 




Wl /n 

+0 [( fc 2 ) 3 ] 


dw i/n 


dfc 2 


k 2 

k^O 


1 d 2 cu 2 


i/ii 


2 d (fc 2 


(k 2 y 


k—>0 


( 10 ) 


so these two dispersion relations are in the desired form of 
Eq 0. Note wjy n = w 2 ^ II (fc 2 ) only permits even powers 


Beyond the hydrodynamic limit, but imposing G1 and 
G2, the equations CD and 0 simplify to: 


d^f/n 


d k' 2 


A /11 


1 ± 1 


— c n 


tr[AG 2 ], 


(15) 


d {k 2 


k^O 


k^O 


tr[Z 2 ]=btr[Z 2 ] 


±tr[G 2 ] — tr[X 2 ] + 


tr[AG 2 

tr[n 2 } 


(16) 


To enforce G1 and G2 the effective coupling between 
the hyperfine states has to vanish, Uab = 0. This 
can be done by imposing a particular transition rate 
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A = — 2^pao Pbo Uab- In addition, fix the hydrody¬ 
namic speed of sound to be 


c 


2 _ 
o — 


ms pao Uaa + mAPBoUBB + Uab {pap^a + pboitt-b) 

2 toaTOb 


(17) 

While one eigenfrequency always remains zero, o/o,i = 0, 
for the second phonon mode we get 


2 4U A B(PAom B +PBom A )cl . . 

u 0,11 ~ 72 ' ( 1 ®) 

n 

Thus the AM corresponds to one massless particle m\ = 0 
and one massive particle mjj = H 2 uJq ii /cq, propagating 
in the acoustic Minkowski space-time in the hydrody¬ 
namic limit. This mass is much smaller that any average 
of the atomic masses if we set Uab Uaa + Ubb- (Al¬ 
though not relevant to the aim of the present paper, we 
stress that such a regime is potentially achievable in an 
experimental setting.) For higher wave numbers we ob¬ 
tain LIV in the form of equation 0 , and the coefficients 
ij -2 and rj 4 for the two modes are: 


%,1/H 

V4,I/II 


f mi/n \ 2 
\M lw ) 

i; 


/ quasiparticle mass \ ‘ 
\ effective Planck scale J 


; (19) 
( 20 ) 


where Mliv = \JrriAmB is defined as the scale of Lorentz 
violations — which is our analogue Planck scale. 

It is quite remarkable that the quadratic coefficients 
m are exactly of the form postulated in several works on 
non-renormalizable EFT with LIV (see e.g. Q). They are 
indeed the squared ratio of the particle mass to the scale 
of Lorentz violation. Moreover we can see from that 
there is no further suppression — after having pulled out 
a factor (h/M^iy ) 2 — for the quartic coefficients 774 , 1 / 11 - 
These coefficients are of order one and generically non- 
universal, (though if desired they can be forced to be 
universal by additional and specific fine tuning). 

Discussion: The suppression of 772 , combined with the 
non-suppression of 774 , is precisely the statement that 
the “naturalness problem” does not arise in the current 
model. We stress this is not a “tree level” result as the 
dispersion relation was computed directly from the fun¬ 
damental Hamiltonian and was not derived via any EFT 
reasoning. Moreover avoidance of the naturalness prob¬ 
lem is not directly related to the tuning of our system 
to reproduce SR in the hydrodynamic limit. In fact our 
conditions for recovering SR at low energies do not a 
priori fix the the 772 coefficient, as its strength after the 
“fine tuning” could still be large (even of order one) if 
the typical mass scale of the massive phonon is not well 
below the atomic mass scale. Instead the smallness of 772 
is directly related to the mass-generating mechanism. 

The key question is now: Why does our model es¬ 
cape the naive predictions of dominant lowest-dimension 
Lorentz violations? (In fact in our model for any p ;§> 


to// the fc 4 LIV term dominates over the order fc 2 one.) 
We here propose a nice interpretation in terms of “emer¬ 
gent symmetry”: Non-zero A simultaneously produces a 
non-zero mass for one of the phonons, and a correspond¬ 
ing non-zero LIV at order fc 2 . (Single BEC systems have 
only fc 4 LIV as described by the Bogoliubov dispersion 
relation.) Let us now drive A —> 0, but keep the condi¬ 
tions Cl and C 2 valid at each stage. (This also requires 
Uab —> 0.) One gets an EFT which at low energies 
describes two non-interacting phonons propagating on a 
common background. (In fact 772 —*• 0 and c/ = c// = Co.) 
This system possesses a SO( 2) symmetry. Non-zero laser 
coupling A softly breaks this SO( 2), the mass degeneracy, 
and low-energy Lorentz invariance. Such soft Lorentz vi¬ 
olation is then characterized (as usual in EFT) by the 
ratio of the scale of the symmetry breaking to//, and 
that of the scale originating the LIV in first place Mliv- 
We stress that the SO( 2) symmetry is an “emergent 
symmetry” as it is not preserved beyond the hydrody¬ 
namic limit: the 774 coefficients are in general different 
if ttia 7 ^ tob, so SO(2) is generically broken at high en¬ 
ergies. Nevertheless this is enough for the protection of 
the lowest- order LIV operators. The lesson to be drawn 
is that emergent symmetries are sufficient to minimize 
the amount of Lorentz violation in the lowest-dimension 
operators of the EFT. 
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